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Gamma-ray binaries are orbital modulated gamma-ray sources in the Galaxy detected both at GeV and TeV 
energies. The high-energy radiation may come from the interaction of energetic electrons injected by a young 
pulsar and photons from the massive companion star. We present a model for the production of high-energy 
radiation where emission, absorption and pair cascading are considered. New observations of LS I +61°303 
by the Fermi Space Telescope revealed an exponential cut-off in the spectrum at a few GeV, inconsistent 
with gamma-gamma absorption. Electrons radiating at GeV and TeV have probably two different origins. We 
investigate whether the emission from the unshocked pulsar wind explains the GeV component in LS I +61°303. 



1. Introduction 



Three binary systems have been discovered to pro- 
duce non-thermal radiation up to very-high energies 
in the Galaxy (see also the contribution of J. Holder 
in these proceedings) namely 1 : LS I +61°303 0,0], 
LS 5039 [j and PSR B1259-63 [5]. HESS J0632+057, 
an unidentified gamma-ray source serendipitously dis- 
covered by HESS may be the fourth system known 
These systems are composed of a Be or O compan- 
ion star and a compact object in an eccentric orbit. 
Except for PSR B1259— 63 where a 48 ms radio pulsar 
has been detected [8J, the nature of the compact star 
remains unknown in the other two binaries. 

These systems are known today as "gamma-ray bi- 
naries" as their non stellar luminosity peaks in gamma 
rays @ . The main feature of these objects is proba- 
bly the orbital modulation of the radiated gamma-ray 
flux. This property has been first uncovered at TeV 
energies by Cherenkov telescopes [T3, El, EH- The 
Fermi Space Telescope has recently reported the first 
detections of the orbital modulation at GeV energies 

EE 

03 (see Richard Dubois' proceeding in this vol- 
ume for more details). 

In this proceeding, we report on our modeling inves- 
tigations of the high and very-high energy radiation 
produced in gamma-ray binaries in the framework of 
the pulsar wind nebula scenario. In this model, the 
non-thermal radiation is emitted by a cooling plasma 
of relativistic electron-positron pairs injected and ac- 
celerated by a young pulsar (see Section 2). We 
present in Section 3 our modeling of the GeV and 
TeV modulation in LS 5039 where emission, absorp- 
tion and pair cascading are considered. In Section 4, 
we discuss the new observations of LS I +61°303 by 
Fermi and investigate whether the emission from an 



1 Note that the detection of Cyg X-3 by Fermi has been 
anounced at the Symposium (see S. Corbel's proceeding). Note 
also that TeV emission from Cyg X-l has been reported by 
MAGIC Q. 



unshocked pulsar wind can explain the GeV compo- 
nent. 



2. The pulsar wind nebula scenario: the 
big picture 



Maraschi & Treves [15[ suggested that the radio to 
gamma-ray radiation could be powered by a young 
(~ 10 4 - 10 5 yr) fast-rotating pulsar in LS I +61°303 
(scenario originally applied to Cyg X-3 [Hj]). The 
same idea has been recently revisited and proposed 
to be at work in the other gamma-ray binaries [9j . In 
this scenario (see the sketch in Fig. [T]), the pulsar in- 
jects energetic electrons-positrons pairs in a relativis- 
tic wind. Particles propagate up to the termination 
shock where the pulsar wind momentum is balanced 
by the ram pressure of the massive star wind. The 
stellar wind can possibly confine the pulsar wind and 
a comet-like structure develops, spiraling around the 
system due to the orbital motion. At the shock, lep- 
tons are randomized, re-accelerated and produce non- 
thermal radiation. 

In the 'shocked' wind, pairs emit synchrotron radia- 
tion and up-scatter the soft photons from the massive 
star to high energies by inverse Compton scattering. 
In the 'unshocked' wind region, pairs radiate only via 
inverse Compton scattering since upstream the termi- 
nation shock, the magnetic field is frozen into the flow 
of particles. No synchrotron radiation is produced. 
The annihilation of gamma rays with the stellar pho- 
tons can also be of major importance in some tight 
binaries such as LS 5039 (see next section). The pro- 
duction of gamma rays by inverse Compton scattering 
is not isotropic in gamma-ray binaries because the soft 
radiation field set by the massive star is anisotropic 
with respect to the location of electrons. Hence, the 
gamma-ray flux depends on the relative position of 
both stars and the observer (Fig. |T|) . 
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Figure 1: This diagram depicts the pulsar wind nebula scenario in gamma-ray binaries. The massive star wind collides 
and confines the relativistic pulsar wind in a collimated outflow. Non-thermal radiation originates from the 'unshocked' 
and 'shocked' part of the pulsar wind. The observed gamma-ray flux depends on the position of both stars with respect 
to the observer (viewing angle tp). 



3. Modeling the high-energy radiation in 
LS 5039 

LS 5039 is composed of a 06.5V star in a compact 
3.9 day orbit with its compact object. The orbital pa- 
rameters are known with good accuracy (the latest so- 
lution can be found in [ijj )• The absence of a Be equa- 
torial wind component makes this system ideal for 
modeling. TeV observations of this binary by HESS 
have shown a stable modulation of the gamma-ray flux 
with a maximum close to inferior conjunction and a 
minimum at superior conjunction [l(J. The main be- 
havior of this modulation can be accounted for by the 
effect of photon-photon annihilation 7 + 7 — > e + + eT 
(see e.g. [18]). In the pulsar wind nebula scenario, the 
high-energy radiation has two distinct contributions. 

The first component originates from the shocked 
wind where pairs are assumed to be isotropized and 
injected with a power-law energy distribution cooling 
down via synchrotron radiation and inverse Compton 
scattering. Combining emission and 77-absorption of 
gamma rays, HESS observations can be reasonably 



well explained. Some physical quantities at the shock 
can be tightly constrained such as the magnetic field 
B = 0.8 ±0.2 G, the spectral index of the injected par- 
ticle distribution p = 2 ±0.3 and the pulsar spin down 
power Lp = 10 36 erg.s _1 (see [l9| for more details). 

This simple model can reproduce most of the spec- 
tral and temporal features at high energies but un- 
derestimates the flux observed by HESS at orbital 
phases where the gamma-ray flux is highly absorbed, 
i.e. close to superior conjunction. The development of 
a cascade of p airs has been proposed to solve this dis- 
crepancy [lOj • The large density of pairs produced by 
photon-photon absorption in the system can reprocess 
a significant amount of the absorbed power. If the am- 
bient magnetic field is high enough, pairs are locally 
isotropized and radiate in all directions, the cascade is 
three-dimensional [20] . But the magnetic field should 
not exceed a few Gauss (5 G in LS 5039) otherwise 
particles will emit preferentially synchrotron radiation 
whose energy is too small to create new pairs in the 
system. Hence the cascade is quenched. 

We present in Figs. the computed lightcurves 
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inates from the unshocked pulsar wind. The inverse 
Compton scattering on a mono-energetic radial wind 
of relativistic pairs produces a sharp line-like compo- 
nent centered at an energy set by the injected Lorentz 
factor of the wind (see [2l[ for more details). In 
LS 5039, this contribution should have a strong spec- 
tral signature but is not observed at high energy sug- 
gesting that the Lorentz factor of the wind might be 
low (< 10 4 ) or very high (> 10 7 ). Alternatively, the 
pulsar wind might have a 'striped' structure and re- 
main highly magnetized up to the termination shock 
(see for instance [22[)- Another scenario has been pro- 
posed where the very-high energy radiation is domi- 
nated by the unshocked pulsar wind [231 ] . 



Figure 2: Computed lightcurve in the Fermi energy band 
(flux > 100 MeV) for LS 5039. The dotted line indicates 
the orbital phase at superior conjunction. The 
inclination of the orbit is i = 40° . 
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Figure 3: Computed lightcurve in the HESS energy band 
(flux > 1 TeV) for LS 5039. The red line gives the total 
flux with the cascade contribution. The blue line shows 
the absorbed primary flux for comparison. Lightcurves 
are averaged over an orbital phase bin A<f> = 0.1 and are 
compared with HESS data points taken from [jc| . The 
dotted line indicates the orbital phase at superior 
conjunction. The inclination of the orbit is i = 40°. 



in the Fermi and HESS energy bands where the full 
calculation of a three-dimensional pair cascade radia- 
tion is taken into account using a Monte Carlo code 
(Cerutti et al., in preparation). The cascade contribu- 
tion is significant at every orbital phase and dominates 
close to superior conjunction, giving a consistent flux 
with HESS observations (Fig. [3]) . The anti-correlation 
between the GeV and the TeV lightcurves is due to 
pair production (see for more details). 

The second high-energy radiation component orig- 



4. New modeling challenges: 
LS I +61°303 seen by Fermi 

The Fermi collaboration has recently reported the 
detection of the gamma-ray binary LS I +61°303 and 
the first detection of an orbital modulation at GeV 
[l3|. The lightcurve peaks close to periastron and is 
minimum close to apastron. The measured spectrum 
is a power- law with a 2.2 photon index and an ex- 
ponential cut-off at 6.3 GeV. This energy cut-off is 
too low to be due to pair production on stellar pho- 
tons. Electrons radiating at GeV and TeV have prob- 
ably two different origins. It is worthwhile to note 
that this spectrum is very similar to those observed 
in gamma-ray pulsars. Fermi might be observing the 
magnetospheric emission from the pulsar in the bi- 
nary system. In this case, the origin of the orbital 
modulation remains unclear. 

Alternatively, the high-energy radiation component 
might come from the unshocked pulsar wind (for an- 
other alternative solution, see the contribution by M. 
Chernyakova in these proceedings). We investigate 
here whether this possibility could provide a viable ex- 
planation for the Fermi observations. Figure [4] shows 
the computed inverse Compton radiation produced by 
the interaction of the stellar photons with the pairs in 
the unshocked wind. Pairs are injected with a soft 
power-law distribution with an index « 3.1 and an 
exponential cut-off at ~ 25 GeV. To account for the 
observed flux, the pulsar luminosity should be high 
L p > 10 37 erg.s" 1 . 

Figure [5] gives the expected GeV modulation using 
the orbital parameters found in [1 71 ] . The lightcurve 
peaks just before periastron, at an orbital phase (</> ~ 
0.15) where both the soft photon density and the an- 
gle of interaction between the pairs and the photons 
is high (close to head-on collision). The peak is fol- 
lowed by a steep decline of the gamma-ray flux. The 
lightcurve is minimum just after periastron (<f) ~ 0-3) 
since the interaction between the pairs and the stel- 
lar photons occurs almost rear-end, hence highly in- 
efficient even though the soft photon density remains 
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Figure 4: Orbit-averaged spectral signature of an 
unshocked pulsar wind in LS I +61°303 (blue line). 
Observations are overplotted for comparison namely 
EGRET (grey bowtie, Fermi (black crosses, |l3j]) 

and MAGIC (red bowtie, 0). 



5. Conclusion 

Leptonic model can explain the general features ob- 
served at high and very-high energy in gamma-ray bi- 
naries. In LS 5039, the GeV and TeV modulation 
are well reproduced but the interplay between emis- 
sion, absorption and pair cascading is complex. Fermi 
measurements of the spectrum and the modulation in 
LS I +61°303 challenge simple models. Particles radi- 
ating at GeV and TeV energies may have two different 
origins. The GeV component might be the signature 
of an unshocked pulsar wind in the system. More in- 
vestigations are required to understand the origin of 
the GeV modulation in LS I +61°303. 
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Figure 5: Theoretical lightcurve (flux > 100 MeV) 
computed in LS I +61°303 (blue solid line). The 
gamma-ray flux is assumed to originate from the 
unshocked pulsar wind. The dotted lines show the orbital 
phase at periastron (blue) and apastron (red). Orbital 
parameters are taken from [l?]]. The folded lightcurve 
measured by Fermi [Tj|] is overplotted for comparison 
(black solid line). 



high at this phase. Nonetheless, this model cannot 
account for the observed GeV modulation. A signif- 
icant shift in phase (A</> ss +0.25) of the theoretical 
lightcurve would be needed to match observations, but 
there is no obvious physical motivation for this yet. 
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